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We study the fate of the energy deposited by a jet in a heavy ion 
colhsion assuming that the médium created is opaque (jets quickly lose en- 
ergy) and its viscosity is so low that the energy lost by the jet is quickly 
thermalized. The expectation is that under these conditions the energy 
deposited gives rise to a Mach cone. We argüe that, in general, the be- 
havior of the system is different from the naive expectation and it depends 
strongly on the assumptions made about the energy and momentum de- 
posited by the jet into the médium. We compare our phenomenological 
hydrodynamic calculations performed in a static médium for a variety of 
energy-momentum sources (including a pQCD-based calculation) with the 
exact strong coupling limit obtained within the AdS / CFT correspondence. 
We also discuss the observability of hydrodynamical features triggered by 
jets in experimentally measured two-particle correlations at RHIC. 
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1. Introduction 



One of the most prominent experimental discoveries made at the Rel- 
ativistic Heavy Ion Collider (RHIC) has been the suppression of highly 
energetic partióles [H El (Sj |4| , which suggests that the matter created at 
RHIC is a color-opaque, high density médium of colored degrees of freedom 
where fast partons quickly lose energy by gluon emission ¡5l [H [71 [U [9l [lÜl [TTl 
[T2I [T3l [H] . Measurements of anisotropies in soft particle momentum distri- 
butions [H El [3111] have further indicated that soft degrees of freedom are 
approximately thermalized. The degree of thermalization has been found 
to be considerably above the predictions [15j obtained within perturbative 
Quantum Chromodynamics (pQCD) and, in fact, it seems to be compatible 
with the "perfect fluid" scenario where the strongly coupled Quark-Gluon 
Plasma (sQGP) has almost zero viscosity [l6l|l71ligil9l|20lEIlE2l[23l[2l¡. 
These two findings, taken together, suggest that the energy deposited by the 
jet into the médium is thermalized and becomes part of the fluid. The the- 
ory governing the further evolution of this energy is that of hydrodynamic 
sound waves. 

It was noticed in the 19th century that, when a probé travels through a 
fluid with a speed greater than the speed of sound, the energy deposited by 
the probé creates a forward moving conical shock-wave. The reason for this 
can be seen in Fig. [1) at each point the deposited energy becomes a sound- 
wave, which moves at the speed of sound Cg. The spherical waves interfere 
coherently creating the shock wave [25j. Fig. [T] can be used, together with a 
simple geometrical argument, to show that the cone angle is related to the 
speed of sound Cg and the velocity of the jet v as 



and, thus, an experimental observation of a Mach cone provides both a proof 
of the ñuid-like behavior of the underlying system and a direct probé of its 
equation of state. This idea actually predates the current investigations 
about sQGP and even quark-gluon plasma by decades }26l I27j . 

Similar ly to ordinary sound waves, Mach cone shocks dissipate expo- 
nentially (~ with respect to the wave-number k and the distance 

traveled x, where the characteristic sound wave attenuation length F is re- 
lated to the shear viscosity ry, the energy density e, and the pressure p 



COSCpM = — 
V 



(1) 



F = - 



4 r] 



(2) 



3 e 4-p 



Therefore, the presence of a Mach cone signal would be an exponentially 
precise confirmation of the low viscosity fluid limit, which is in principie 
more sensitive to rj than even anisotropic flow. 
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As with all heavy ion observables sensitive to hydrodynamic behavior, 
the Mach cone signal suffers from the problem that we do not "see" the 
ñuid directly: only the final many-particle correlations are measured and 
they are sensitive to all stages of the hydrodynamic evolution including 
the late (presumably non-thermalized) stages and freeze-out (which is not 
understood from first principies). A rough approximation is to assume that 
at a certain locus in space-time S'^ = (í, x) the mean free path goes from 
zero to infinity. This locus can be defined in terms of a local criterion 
(e.g. a common freeze-out temperature) , or using a simple global geometry 
(isochronous freeze-out). Using Stoke's theorem, as well as entropy, energy 
and momentum conservation yields the famous Cooper-Frye (CF) formula 
[28] 

E^ = j d^^^P^^f{U^P^,T) (3) 



4 ZAKOPANEV4.1 PRINTED ON JANUARY 6, 2009 




O 71 (|) 

Fig. 2. (Color online) A schematic representation of a Mach conc in hcavy ion 
collisions. 

where is the collective flow vector, f{E, T) the standard Boltzmann, 
Bose-Einstein or Fermi-Dirac distribution function (in terms of tempera- 
ture T), and is the 4-momentum vector {E,p) oí the associated particle. 
Using an azimuthal coordínate system and putting the jet direction at the 
origin yields a characteristic distribution dN/dcp with two Mach cone-hke 
peaks shown in Fig. [2j Note that "theoretically" this is an average distri- 
bution since we define the near-side jet to he at (j) = 0. 

Experimentally, however, this is a 2-particle correlation: the experiment 
measures a high momentum "jet" particle (the near-side trigger) and then 
looks at the correlation between the trigger and softer (sensitive to ñow) 
particles in the opposite direction (away-side región, where it is assumed 
that the jet passed through the médium, was suppressed, and its energy 
thermalized into Mach shocks). Tantalizingly, something similar does seem 
to be observed in the experiment (Fig. [3]): when a hard (jet-like) particle 
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Fig. 3. (Color online) Two-particle correlations induced by jets in d-Au coUisions 
(red) and Au-Au coUisions (blue) at RHIC [29l [301 El] ■ With a hard away-side 
trigger (left, p?j?''°'^ > 2 GeV/c), the away-side peak disappears in agreement with 
the hypothesis that the jet is absorbed by the médium. However, if the away- 
side trigger momentum is lowered (right panel, 0.15 > p?j?^°'^ > 4 GeV/c), the 
peak reappears (as expected from momentum conservation) and shows a cone- 
like pattern. The absence of similar correlations in the d-Au "control" experiment 
confirms that this is an efFect of the "médium" rather than a deviation of the initial 
conditions. 

is correlated with a soft (medium-like) particle, the "missing" away-side 
jet signal reappears and shows a structure very similar to that of a Mach 
cone with an angle cióse to the expectation for the ideal gas speed of sound 
Cs = l/\/3. 

It should be said that this result has generated a lot of skepticism and 
some debate (for a review see [291 EOl El]). The "dip" between the Mach 
cone peaks arises most clearly when the background 2-particle correlation 
(arising, for example, from the ellipticity of the initial fireball) is subtracted 
(ZYAM method). This subtraction makes possibly unjustified assumptions, 
e.g., that the Mach cone correlation and the elliptic correlation are inde- 
pendent. Moreover, different theoretical interpretations have been given to 
the apparent conical structure [521 [331 IMl [3S] > which indicates that further 
studies involving, for instance, 3-particIe correlations p6] and the depen- 
dence of the Mach cone angle on various other variables [ST] [38] should be 
pursued. 

Nevertheless, the experimental observation of something that looks like a 
Mach cone has greatly enhanced the theoretical interest on this phenomenon 
[Ml[10l[3Illl2l[l3l[l5llia[l6l[13[l8l[l9l[50l|5^ 

[59] [ÜOl W\\ [62] . In particular, an interesting question that still remains 
to be properly answered is to what extent does the "naive" Mach cone 
picture introduced earlier survive when realistic physics conditions are in- 
troduced. Recent studies seem to indícate the Cone-|-diffusion wake pic- 
ture, and the relative weight of the two, are remarkably independent of the 
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energy-deposition scenario 

In the next few sections we shall examine this issue in detail. Sec- 
tion [2] describes an analysis of Mach cone propagation using ideal three- 
dimensional full (non-linear) hydrodynamics, which shows that the observ- 
able signal depends crucially on how energy and momentum are deposited 
by the jet into the médium. Section [3] uses the Anti-de Sitter/Conformal 
Field Theory (AdS/CFT) correspondence to try to obtain insights into the 
non-equilibrium dynamics displayed by this problem, introduces the "Head 
+ Neck" decomposition of the space-time región near the jet (based on the 
degree of thermalization of the médium), and derives the contribution of 
each región to the partióle correlations. Section U] does an analogous anal- 
ysis with an energy-momentum source calculated in pQCD coupled with 
ideal hydrodynamics and describes the observable differences with respect 
to results obtained using AdS/CFT. 

We use natural units and the Minkowski metric ^^¡y = diag(— , -|-, + , +). 
Lorentz Índices are denoted with Greek letters = O, . . . , 3. Throughout 
this paper we shall use a coordínate system corresponding to axial symmetry 
with respect to the trigger jet axis. The components of a generic vector x 
will be 

/ t \ 

Xl = X — vt 
X± eos if 

y x± sin Lf j 

coordinates, where are taken to be parallel/perpendicular to the near- 
side jet. The jet is assumed to be moving transversely with respect to the 
beam (z axis) so the rapidity of the jet is taken to be zero. The away-side 
direction in this system is (0, 1, 0, 0). 



(4) 



2. Mach Cones in Hydrodynamics 

The geometrical argument for Mach cone formation and the angle for- 
mula in Eq. ([1]) are only valid in the linearized hydrodynamics limit where 
the energy of the sound waves is small compared to the energy density of 
the background. The physical applicability of this condition is doubtful 
since the order of magnitude of the jet's "size" ~ 1 fm, which leads to an 
overwhelming energy density cióse to it even for modérate rates of energy 
deposition. A "pile-up" of sound waves in front of the jet as it travels 
through the médium can also result in strong non-linear corrections to the 
ñuid's ñow profile. Furthermore, if the energy deposition is small with re- 
spect to the background, an observable signal will not survive the thermal 
background ñuctuations inherent in a CF freeze-out. It is relatively easy 
to show analytically ¡Hl [56] that this is the case unless either the strength 
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Fig. 4. (Color online) Temperature (left panel) and ñow (riglit panel) profiles for a 
hydrodynamic simulation of a jet with difFerent energy and momentum deposition 
schemes in a static médium (from Ref. [53]). 



of the signal is large enough to require a non-linearized treatment or the 
momentum is large enough to put the thermalization assumption in doubt 
ñ 

For associated (massless) partióles with = {pt,Pt cos{'K—(j)),pT sin(7r— 
(/)),0) the momentum distribution at mid rapidity is given by 

= / dS^P'^ [/(C/^ T) - (PO, To)] (5) 

^ "Harder" partióles are more contaminated by non-thermal processes such as minijet 
fragmentation. While this is "obvious" since the effective mean free path grows with 
the partióle momentum, the scale at which the médium stops being the main source 
of partióles is not known preoisely (though assumed to be in the pt ~ 1 — 2 GeV 
range). 
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We subtract a non-ñowing = (1,0) thermal T = Tq constant back- 
ground yield via feg = /|c/ai=o,t=To- Viscous corrections to the Boltzmann 
distribution function [6lj produce subleading contributions that are negligi- 
ble in the linearized approximation. Choosing an isochronous ansatz where 
dTi^ = x±dx±dxidíp (1,0,0,0), the Boltzmann exponent can be expanded 
up to corrections 0{{U)'^). The associated away side azimuthal distribution 
at mid-rapidity f{pT,<P) = dN /pTdpTdyd(¡)\y=o with respect to the nuclear 
beam axis is then given after integrating over (/? by 

f{PTA)='^'^PT / íÍ2;i(ÍX_LX_LX (6) 

Jt, 

(exp {-^ [Uo - Ui cos(7r - 0)]} lo(ax) - e-P^/^°) 

where a± = ptU± sin(7r — (j))/T and Iq is the modified Bessel function. In 
the linearized approximation a± <^ 1 and, thus, we can use the expansión 
for the Bessel function 

lim/o(x) = l + ^ + 0(x4) (7) 



to get the approximate equation for the distribution 



(8) 



Note that deviations from isotropy are then controlled by the following 
global moments (AT) = ^^dx\dx ^ I^T and (C/i) = ^^dx\dx^x^\J\. 

It is clear that in the strict linearized limit the azimuthal distribution 
only has a trivial broad peak at <?í> = tt. A double-peaked structure in the 
away-side of the jet correlation function can only arise when the Bessel 
function expansión is invalid, i.e., away from the linearized limit or for large 
Pt ^ Tg. For such large momenta, contamination from non-thermalized 
degrees of freedom and coalescence [651 ESI EZl [MI [ÜSl [7D1 |7l] effects are 
non-negligible. 

On the other hand, for a non-linear perturbation in the médium it is not 
clear that the Mach cone angle is anything like the one derived geometrically 
in the Introduction. For large energy depositions, the more appropriate de- 
scription is that of an angular shock, i.e., a step function in energy density. 
This problem was analyzed analytically in ¡72j . where the angle was found 
to be in general larger than the geometrical expectation from Mach's law. 
Deviations from equilibration, expected in the región where the energy de- 
position from the jet is comparable to the local energy density, may also 
lead to a different cone angle. 
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In addition, the Mach cone is not the only collective mode that can be 
excited by the moving jet. The ratio of momentum to energy deposition is 
not known but a common assumption is to have the same amount of en- 
ergy and momentum being deposited by the jet. Longitudinal momentum 
deposition results in a column of ñuid that flows in the direction of the 
jet's motion. This structure is known as a the "diffusion wake". We have 
taken ah these effects into account in Ref. ¡53j where we studied energy and 
momentum deposition in a static médium using the 3D ideal hydrodynam- 
ics SHASTA code presented in [73]. Note that the simulations presented 
below only refer to a static thermal background, and as such neglect the 
effects of transverse and longitudinal ñow (which are a simple deformation 
in the linearized approximation \42\ but can be non-trivial in full nonlinear 
hydrodynamics [63j). 

Fig. [3] shows what happens when the momentum deposition is included 
in the full hydrodynamic simulations [53]: the temperature profile remains 
invariant and maintains the correct angle. However, the ñow profile acquires 
an additional component co-moving with the jet. See Fig. [5] for a sketch 
of all these effects combined. It then becomes clear that the appearance 
of a cone-like signal in hydrodynamics is not assured, since the strongest 
signáis may come from contributions which are not conical and possibly 
not locally equilibrated. More complicated momentum depositions, such as 
including transverse momentum, or introducing a momentum dependence 
of the deposited energy (the so-called "Bragg peak") do not change these 
conclusions qualitatively ¡63] . 

-Fig. [ü] shows what happens after freeze out: when conical distributions 
domínate, the system does exhibit a cone signal with the correct angle, al- 
beit the momenta of the associated particles are so high that one is likely to 
question their effective thermalization once the ideal fiuid approximation is 
relaxed. On the other hand, when the diffusion wake is significant the only 
correlation visible is a unique peak in the away-side, which is indistinguish- 
able from a generic peak expected from momentum conservation. 

To strengthen this conclusión, note that the puré energy deposition 
scheme may be completely undetectable in a realistic and approximately 
isothermal freeze-out condition since hotter regions will simply freeze-out 
later. Momentum ñow, however, will persist independently of freeze-out 
and will in general be modified only a little by the last (cooler) stages of the 
fireball evolution. Thus, it is the right panel of Fig. [J] that shows the ñuid 
correlations most apt to be imprinted on two-particle correlation functions. 

In sum, a phenomenological hydrodynamical approach shows that Mach 
cones are not guaranteed to appear in the final angular correlation functions. 
The problem of the away-side jet correlations in the sQGP is a very Inter- 
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Fig. 5. (Color online) A schematic representation of a Mach cone solution including 
non-linear and non-equilibrium effects. 



esting and complicated subject because of its inherent out of equilibrium, 
non-linear, and non-perturbative features. An interesting development in 
this direction was provided by the possibility of using string theory dualities 
to calcúlate observables in a strongly coupled gauge theory without the "a 
priori" equilibration assumptions. 

3. Mach Cones in AdS/CFT 

3.1. Introduction 

The energy-momentum tensor of a system composed of a heavy quark 
passing through a strongly coupled A » 1 (A = Qsym^c is the t'Hooft 
coupling) = 4 SU{Nc) Yang-Mills plasma at finite temperature T can 
be computed using the AdS/CFT correspondence [36 l liTl HHl H9 l [50 l [51 ] . 
first conjectured in [Tlj. According to the correspondence, gauge invariant 
observable quantities in a strongly coupled = 4 SYM theory can be de- 
termined using weakly coupled 10 dimensional type IIB superstring theory, 
where 5 of the dimensions are Anti-de Sitter and the other 5 correspond to 
a 5-sphere. 
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Fig. 6. (Color onlinc) EfFect of freeze-out for various energy-momentuni depositions 
[55] . Note that angle (¡> is shifted in the left-hand plot in such a way that the away- 
side peak is located at </> = 0. 



The problem of a heavy quark moving at constant speed in a strongly 
coupled finite temperature = 4 SYM médium can be analyzed (Fig. 
[7|) by considering metric ñuctuations due to a string that is hanging down 
from the boundary of an AdS Schwarzschild (AdS-SS) background geometry 
[75l [Tül [77| . Quantum ñuctuations can be neglected for a slowly moving 
heavy quark [75]. In this hmit, the action that describes the supergravity 
approximation to type IIB string theory in an AdS-SS background and the 
classical string is given by the sum of the fohowing partial actions 

and 

Ang = / \/-Gf}daXt^dpX-éa, (10) 



^ Note that "Mach-like" signáis found in experiment so far are triggered by light 
quark/gluon jets P l4l [29l |3Ü] . 
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Fig. 7. (Color online) The Mach cone set-up associated with heavy quarks within 
the AdS/CFT correspondence (the figure on the left-hand side was taken from 
Ref. ¡46¡). The finite-temperature médium corresponds to an event horizon and 
the quark to the tip of a fundamental string stretched between the horizon and the 
boundary of AdS^. 



where L is the radius of AdS^, G5 = ttL'^/2N^, a' = L^/\/A, 



total metric, and G^/^u is the metric of the unperturbed AdS-SS black hole 
(without effects from the string), which can be obtained from 



G^u is the 



ds' 



z2 



-g{z)df + df^ + 



dz' 
W) 



where z goes from O at the AdS boundary to zq = l/vrT at the black hole 
horizon (T is the Hawking temperature associated with the black hole) and 
g{z) = 1 — (z/zq)^- The string coordinates X^^{a,T) in the Nambu-Goto 
action in Eq. (¡lOp are chosen in such a way that the string endpoint (which 
corresponds to the heavy quark in the 4-dimensionaI boundary) moves at 
constant speed v and no energy ñows from the horizon into the string ¡75[ 

Eü]. 

Minimizing the action S with respect to the metric G leads to the full 
set of Einstein's equations. It is sufñcient for our purposes here to consider 
instead the linearized Einstein's equations for the metric ñuctuations 

.(0) 



which are defined via G^y = G]^J + h^y. It can be shown ¡3ül [1T| that th_e 



contribution from the moving quark to the total energy-momentum tensoi0 
is Tquark = ■^\J Qi where the tensor Q is obtained by expanding h in 
powers of z near the boundary, i.e., h ^ Q z^. An example of the formidable 



analytical power of AdS / CFT calculations was given by Yarom in Ref. 
(see also ¡78] ) where he computed the total energy-momentum tensor in the 



By Tf¿u we mean {Ti_i^), though we will drop the (...) notation for brevity. 
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lab frame that describes the near-quark región 

Tj, = Po diag{3, 1, 1, r2} + ^ Pq AT^,(xi, r) (12) 
where the explicit form of AT^i, is 

. v{r^{-b+lZv^-Sv*)+{-b+llv^)xl)xi , . 

Ai+í = a— 7^7:5 — , lú 



' 24(r2(l-^2)+x2)''^" ^ ' 

_ (1-1.2 )t,2(llxf+8r2(l-^2))r 

¿Xltr - -a — — — 5/2 , (i5j 



i>(r2(8-13D2+5»;*)+{ll-5»;2)x2)xi , . 

^^^1^1 70/ 2n 2AV2 (16) 

72[r^{l—v^)+xf] 



^ 72(r2(l-t,2)+a;2)"/2 ' ^ 



l,(l-i;2)(5^2(i_^2)+8x2)xi 

ATrr = -a TTsTr^ — , (18) 

72{r^{l-v^)+xlf' 

and a = ^^Txr'^ . In the equations above we used dimensionless coor- 
dinates in which the distance is normahzed to 1/ (vrT) and also defined 
^ = and 7 as the quark time dilation factor l/Vl — . More- 

over, Po = iVcVr^T^/S + ©(TV^) is the pressure of the ideal SYM plasma 

Note that it is assumed throughout the derivation of Eq. (fT2]) that the 
metric disturbances caused by the moving string are small in comparison 
to the AdS^ background metric. Therefore, this result is correct as long 
as this condition is fulfilled. In fact, since AT^i, scales inversely with the 
total distance from the quark, the región where the condition ^Ar^,^ < 1 
(or, equivalently, h small in comparison to GW) holds can be taken to 
be arbitrarily small as long as the limit where Nc ^ 00 and A is large is 
employed. However, in order to evalúate the relevance of this approach to 
heavy ion collisions, one could set Nc = 3, X = 3tt (a^ = 0.25), and 7 = 10, 
which gives ^ > 30. This valué of ^ then sets a lower bound of S/vrT on 
the mínimum distance from the quark that marginally fulfills the condition 
^AT^, < 1. 

We stress that the energy-momentum tensor shown above is not a so- 
lution of the hydrodynamic equations but rather the full non-equilibriuni 
result in the strong coupling limit. The resemblance of this set-up to the 
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hydrodynamic calculations sections [D and [2] indicates that, for strongly cou- 
pled field theories, jet energy deposition is really reduced to something that 
looks hydrodynamical and linear reasonably cióse ~ 5/(7rT) from the jet. 



3.2. Comparison to Hydrodynamics 

The disturbances in the fluid caused by the moving jet are expected to 
behave hydrodynamically in the región sufíiciently far from where the jet is 
presently located. However, in the near zone cióse to the heavy quark hy- 
drodynamics must break down and this can be checked explicitly by looking 
at the isotropy in the Landau frame (denoted henceforward by brackets and 
subscript L, (.-Ol) at each point [56], which is defined by the condition that 
(T(^)l = 0. The boost with respect to the lab frame defines the hydrody- 
namic flow vector U^. Note that, unless the system is a coherent field where 
the phase velocity is equal to the speed of light (such would be the case of 
an electromagnetic wave), this transformation is always possible. This can 
be accomplished by solving a system of two equations for the two space-like 
components of U^, Ui and Ur {Ug = 0), 

(ro^)^ = AfrJ,AÍÍ = o, (20) 

where A^' is a general coordínate dependent Lorentz transformation 

T \ 



and 7 = [/" = \/l + Ü^. Using the representation 



r^ = ( 1 (22) 



where e = TJq and fjj = {i, j = 1, 2, 3), we obtain that Eq. ([2Ü]) becomes 



U 



1 



^(7-1) 



-fS-fU). (23) 



For finite Nc and A this equation can only be solved numerically. However, 
for very large Nc and large A one can use that, to leading order in 1/Nc, 
7 ~ 1, e ~ 3Po) and f ~ Pq- Using these approximations, one can then 
obtain that 

£/»A. (24) 
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In deriving the equation above we used that ^ ~ 0{^/X/N^). The fact 
that \U\ ~ 0{\^/N^) in the large Nc, A hmit imphes that nonUnear terms 
in the hydrodynamic description are subleading contributions that can be 
neglected. Thus, if the system's dynamics can be described by hydrody- 
namics, these equations have to be hnear for the present string theory setup 
obtained in the large A hmit. However, at finite Nc and A nonhnear 
effects are expected to be relevant. These nonhnear effects can only be 
properly taken into account by incorporating subleading 1 /Nc corrections. 



Knudsen Kn(x) v=0.9 




x.,= z-vt {MkTq) 

Fig. 8. (Color online) The local Knudsen number field for the near zone Yarom 
stress. Note the 3Kn > 1 región that defines the Knudsen Neck zone when v = 0.9. 



In hydrodynamics, the Knudsen number Kn is defined as the ratio be- 
tween the mean free path Imfp and a characteristic spatial dimensión of 
the system q. Hydrodynamics is applicable when Kn = Imfp/q ^ In 
conformal field theories at finite temperature, the only dimensionful param- 
eter is given by the temperature T and, thus, both Imfp and q should be 
proportional to 1/T. However, the mean free path is not a well defined 
quantity in = 4 SYM theories at very strong coupling. Nevertheless, one 
can still define an eíFective Knudsen field in terms of the sound attenuation 
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length r and the Poynting vector S as follows [52] 



s 



(25) 



It is easy to check that in the weak coupling limit this definition reduces 
to the usual one corresponding to the mean free path over the macroscopic 
system size (in this case the scale is defined by the flow gradient) or, in other 
words, the number of times a microscopic degree of freedom interacts while 
traversing a macroscopic región the system. Also note that this definition 
is weh defined in the supergravity hmit. As Fig. [8] shows, when v = 0.9 the 
región where hydrodynamics provides a good approximation of the near- 
quark región approximately coincides with the locus defined by the local 
Knudsen number being Kn^^ > 3, which then corresponds to the onset of 
hydrodynamics (Kn « 1) as a sensible approximation. We will later use 
this condition to further analyze our AdS/CFT results and compare it to 
a pQCD-like solution. A direct comparison of the near-zone Yarom tensor 
and a first-order Navier-Stokes ansatz showed ¡52] that a hydrodynamic 
description of the disturbances caused by the heavy quark is valid down to 
distances of about 1/T from the heavy quark. 



3.3. Observability of the AdS/CFT Mach Cone 

The fact that in the supergravity approximation T^*^ can be described by 
(mostly) linearized hydrodynamics means that once the system breaks up 
into particles, a conical signal in the corresponding angular correlations may 
be washed out by thermal smearing, as discussed in the beginning of Section 
[21 This means that a detectable Mach cone-like signal may also come from 
the región where the linearized approximation is not valid. In fact, we will 
show that the only detectable cone-like signal from the AdS/CFT solution 
comes precisely from the región that is not fully thermalized. 

The T^'^ and T^^ components of the energy momentum tensor (com- 
puted within the supergravity approximation) that describe both the near 
región and the far zone were computed numerically by Gubser, Pufu, and 
Yarom in Ref. [lÜ]. In Fig. [ü] we show the energy density perturbation 
Ae{xi, x±)/esYM computed using the data from Ref. [49J due to a heavy 
quark jet with t> = 0.9inaAA = 4 SYM plasma modeled via the AdS/CFT 
string drag model for Nc = 3, X = Qyn-^c = 5.5. The left panel shows the 
far zone (the numbers in the plot label the contours, in per cent as defined 
on the upper-left córner). The Mach wake zone is above the dashed line, 
cos(f)M = ^/W^v), and the Diffusion zone lies below that line. Normalized 
Poynting (momentum flux) vector flow directions are indicated by arrows. 
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The insert shows the nonequihbrium "Neck" zone (with the Coulomb Head 
subtracted) as defined by the condition that Ae/esYM > 0.3. 
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Fig. 9. (Color online) Too (contour) and Tqí (arrows) for the Mach cone AdS/CFT 
solution (the full numerical computation of 09]). The three regions defined in Eq. 
(|26|) are identified. Note the presence of a strong non-hydrodynamic transverse en- 
ergy flow near the core. The dashed line shows the Mach cone line in the linearized 
approximation. 



In order to understand the various regions mentioned above, we decom- 
pose T^y into terms that domínate In diíFerent spatlal scales (see Fig. [9]): 

TM-(x) = T^^ + ¿Tl';,^, + bTZ.^ + JT^l, (26) 

where the far zone "Mach" part of the stress tensor is defined by KrC^ > 3, 
and coincides with the hydrodynamic description, 

6TMach{xi,x^) = jr^ Qf/'^C/- - i^^^^ + ■^d^'^U''^^ - To^'^l 

X e{l-3Kn) (27) 

and the Neck zone is defined by the región cióse to the heavy quark jet where 
the local Knudsen field is large and even uncertainty bounded equilibration 
rates are too small to maintain local equilibrium. The background stress 
tensor in Tq'^ . As shown in [47 1 l56 j the non-equilibrium zone is characterized 
by a stress of the form 

6TNeck{xi,x±) « e{3Kn{x) - 1) , ^^'^(xi, x^) (28) 
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where Y^'^ is a dimensionless "angular" tensor field. At very small distances 
from the jet, ST^eckixi, x±) reduces to the analytic stress tensor defined in 
Eqs. ^MB- 

Within the Neck zone, there is also an inner "Head" región where the 
stress becomes dominated by the contracted Coulomb self field stress of 
the quark 5T^'^^^^. The Head zone can formally be defined as in Ref. [80] 
by equating the analytic Coulomb energy density ¡l6], sc{xi,Xi_)^ to the 
analytic near zone energy density [47| given by Eq. (I13p . This Coulomb 
head boundary is approximately given by 

„2 1^2 2 ^ __}^ (2xj_ + x\f 

^ (7rro)4 74x2(x^/2 + 72x2)2- ^^^^ 

The Head zone is a Lorentz contracted pancake with longitudinal thick- 
ness Axi^c'vrTo ~ 1/7^/2 effective transverse radius Ax^^c'TtTq ~ 

1/7^/2^ which is in agreement with the general considerations in Ref. ¡80¡ . 
However, as we have shown in Fig. [8l for relativistic jets the Neck has a 
two lobe structure. The lobes are nearly independent of v and the lobe 
región thickness is lS.xi^Kn ~ ^/t^Tq » Axi^c". The second thin pancake 
component of the Kn that develops for large 7 (not shown) is similar to the 
shape of the Head zone. The relative independence of the two lobe compo- 
nent of the Neck zone on v is in agreement with the parametric dependence 
Axi^AT oc 2/7rTo ~ GT expected from the bound of dissipation rates imposed 
by the uncertainty principie. 

We now turn to the observable consequences of the AdS string drag 
stress model by assuming the Cooper-Frye hadronization scheme [28j with 
isochronous freeze-out discussed in section [2]0. Our system of coordinates 
is explained in more detall in Fig. [TÜJ As the discussion in section [5] has 
shown, the only way that a nontrivial angular correlation can arise in the 
soft degrees of freedom within the AdS/CFT string drag model is if we relax 
the linearized approximation, i.e., the formal Nc, A ^ 00 but ^/X/N"^ — > O 
condition used to derive the stress and boldly extrapólate towards more 
"physical" parameters to make contact with our QCD world. 

We computed /(</); V) with iV^ = 3, A = 5.5 for v = 0.58, 0.75, 0.90, in a 
static uniform background from the tables of T^^ and T^^ (used in |49j ) . Our 
total CF volume is defined by -14 < Xi (vrTo) < 1, O < (vrTo) < 14, and 
(/? G [0, 27r]. Here we define the head of the jet as the volume where ^ > 0.3, 
which roughly corresponds to the región between — 1 < Xi (vrTo) < 1 a-nd 



* This is a strong model assumption on top of the AdS calculus and wiU need 
much closer scrutiny in the future. We have, however, experimented with differ- 
ent hadronization conditions and observables in [58] and found the results of this 
section, as well as the analogous resuhs in section \A72\ to be quahtatively the same. 
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O < X±_ (vrTo) < 2. We show results for the azimuthal angular correlations 
in Fig. [TTJ The blue curves exelude the chromo-viscous Neck zone from the 
CF volume , the red Neck curves only include the Neck zone approximated 
here by 5T^^{x) > Q.2,esYM- 

As can be seen, only the red Neck curves display the double-peak struc- 
ture while the "Mach" zone is too weak even in the Nc = 2> extrapolation 
of AdS to produce a dip at = vr. For v = 0.9 the two peaks from the 
Neck zone appear at angles accidentally similar to the putative Mach cone 
angle. Particles coming out of the "Neck" región will however probably not 
be thermalized (exponential, with the slope parameter given by tempera- 
ture and ñow as in Eq. [5]), though it is difñcult to say whether one could 
establish "lack of thermalization" conclusively from experimental data. 

The fact that the "cone-like signal" arising from the neck región is not a 
true Mach cone has more direct phenomenological consequences: as shown 
in Fig. [TU the angle between the double peaks is relatively independent of 
the jet velocity v (here tested for v = 0.58,0.75 and 0.99), which violates 
the Mach's law dependence (indicated by the small arrows). We propose 
that looking for deviations from Mach's law for supersonic but not ultra- 
relativistic identified heavy quark jets could test this novel prediction of the 
AdS/CFT drag model. Unlike light quark jets, even high momentum heavy 
quarks move at velocities significantly smaller than c. A sean of Mach cone 




Fig. 10. (Color online) Schematic of the geometry used here. The trigger jet corre- 
sponds to the heavy quark produced near the surface in the —x direction transverse 
to beam axis "z". The away side jet nioves in x direction with velocity v. The 
comoving coordínate is xi — x — vt and the transverse radial coordínate relative 
to this is Xp. A Mach wake (solid blue line) is produced at azimuthal angle ¡pM in 
the xi — Xp plañe. 
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angles with the jet velocity should be, therefore, experimentally feasible. 
These results permit us to use the Mach-cone like signal to understand the 
"phenomenology" of AdS/CFT, since it is not at all obvious that any model 
wiU have a sizeable enough Neck región with the correct ñow behavior. In 
fact, an exphcit counter-example is the pQCD model proposed by ^59j. 



dN((^)/dyd(^ - dN(0)/dyd^ 



0.1 



0.08 



0.06 



0.04 



0.02 



, 1 1 1 , 1 1 1 , 1 


I I I I I I I I I I I 
x30 


1 1 , 1 1 1 , 1 1 1 , 

X 3000 


V = 0.90 


V = 0.75 


v = 0.58 - 






y' " 


' 1 
' 1 


1 1 
1 1 


* * l 
* 1 
* 1 




i' ^, ; 




^/j/ 1 1 1 Va 


M \ / U 

/i t t ,V>L 





123456 123456 123456 



Fig. 11. (Color online) Mid-rapidity azimuthal away-side associated angular dis- 
tribution from the Cooper-Frye freeze-out of the AdS/CFT string drag model 
{T{x),U{x)) fields from ¡49]. Three cases for various heavy quark jet velocity 
and associated hadron transverse momentum ranges, 1 : (v = 0.9, pt/t^To = 4-5, 
2 : {v — 0.75,pt/ttTq — 5-6), and 3 : {v = 0.58,pT/7rTo = 6-7), are compared. The 
short arrows show the expected Mach angles. The red curves showing the double 
shoulder away side dip (conical) correlations are from the Neck región defined here 
as where Ae/esYM > 0.3 (see Fig.[9l). The blue curves result from integrating only 
in the far Mach zone outside the Neck región and show no sign of the weak Mach 
wake seen in Fig.[9]because the NO-GO freeze-out theorem (sec[2]) remains in forcé 
even for our A^c = 3, A = 5.5 downward extrapolation from the supergravity limit. 
The sum total correlation exhibits a double shoulder correlation for v < 0.9 arising 
from the chromo-viscous near zone that is however unrelated to the Mach wakes 
seen in Fig. [Ü) 
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4. Mach Cones in pQCD Coupled to Ideal Hydrodynamics 

4-1. Energy-Momentum Source Computed in pQCD 

As before, we will decompose the energy-momentum tensor into a Head, 
Neck and Mach part according to Eq. (¡26¡). However, in this section we focus 
on the case of a pQCD plasma treated in the chromo-hydrodynamic Umit 
[59] , The background stress tensor in this case is Tq'^ = diag(eo,í'0;í'0)Po)) 
where eo = Stt^Tq/IS is the background energy density of gas of massless 
SU (3) gluons (background temperature Tq) and po = Eq/S. The role played 
by the vacuum contribution to the energy-momentum tensor, 5T^^^^i{X), 
which is associated with the classical, non-Abelian Lorentz boosted Coulomb 
field created by the fast moving parton, is to produce the anomalous re- 
sponse of the médium denoted by the Neck component. Here, we assume, 
as in section [3^ and [80], that the very near bare Coulomb field stress zone, 
in which STj^'^^i ~ 0(l/x^), is the self field stress of the heavy quark and 
does not fragment into associated hadrons. 

The far zone "Mach" part of the stress can be expressed in terms of 
the axially symmetric local temperature T{X) and fluid (Landau) flow ve- 
locity fields U°'{X) through the first-order Navier-Stokes stress form Eq. 
([27|) . Even though we will use the pQCD chromo-viscous source computed 
in linear response in |59i 160], we will assume here the perfect fluid to maxi- 
mize the freeze-out azimuthal conical signature that is otherwise even more 
washed out when viscous dissipation is taken into account. Our aim here 
is not to flt RHIC data but rather to contrast weakly coupled and strongly 
coupled plasma response effects in the most idealized conditions of a uni- 
form static plasma coupled to the external Lorentz contracted color (E", B") 
fields associated with a uniformly moving supersonic color charge. 

Thus, while we assume in the background the perfect fluid r]/s = O limit 
of the full anomalous chromo-viscous equations derived in ¿81j , we also retain 
the anomalous difi'usion stress Neufeld source [59l [60] that we rewrite in the 
more easily recognized Joule heating form (see Eqs. (6.2 - 6.11) of [81]) 

a^T'^^ = cT^ = F^°V;Í = (F^°''C„/3^*F^^'^) (30) 

where F^^"'[X) is the external Yang- Mills field tensor and 

J^{X) = J d^K/{2TT)^exp{-iK ■ X)r{K) (31) 

is the color current that is related via Ohm's law to F^'^"'{K) through the 
(diagonal in color) conductivity rank three tensor C^ap{K). The * denotes 
a convolution over the nonlocal non-static conductive dynamical response 
of the polarizable plasma. 
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The source is then [59] 

= / ^' ^^^^ 

where -Dq/3 is a quadratic form in the external field tensor F^^^"" components 
and /(X, P) is the parton distribution function. In fact, the above and Eq. 
(28) of [59] factorize into the Joule heating (F" • J°-Y form above. 

The covariant generahzation of Neufeld's source is easiest in Fourier 
decomposition with Jy{K) = Cu^a{K)F^°'"'{K) and the color conductivity 
expression derived in [82] 

C,.,{K) = ^9' j d'P pJf¡p'^^^J o{P) (33) 

where fo{P) = {N^ — l) G{P) is the effective plasma equilibrium distribu- 
tion with G{P) = {2-K''^y^e{Po)S{P'^)/{e^'>/'^ -1). Here, is the 4-velocity 
of the plasma as in Eq. (¡271) . For an isotropic plasma C^apiK) = c¡j,a{K)Up. 
In the long wavelength limit, C^°^I^{K ^ 0) = r*m^ g''^"C/^/3, where 
= g^T"^ is the Debye screening mass for a noninteracting plasma of 
massless SU {'i) gluons in thermal equilibrium. 

The relaxation or decoherence time r* is of the generic form noted in 

m 

1111 

- = - + - + — (34) 

"Tp "Te "Tan 

with Tp oc (a^T ln(l/Q5))~^ being the collisional momentum relaxation time 
[151 [83], Te = {oisNcTln{mD/'mM))~^ being the color diffusion time defined 
in [82] in terms of the Debye electric and assumed 0{g^T) magnetic screen- 
ing masses, and Tan oc (mi:)(??| Vf/I/Ts)"*^/^)"^ being the anomalous strong 
electric and magnetic field relaxation time derived in Eq. (6.42) of [81]. Note 
that one can express 

_ 1 1 

^""^ ~ ^Kn{X) ^^^^ 

in terms of the local Knudsen number. However, because r] oc t*sT, Eq. (|35p 
is really an implicit equation for Tan- Combining these relations and taking 
into account the uncertainty principie constraint discussed in section [3l and 
(r*^l/(3r)) m, 

^=T{ai lng-^+ as g^ Ing-^ + 03 g^/Kn) ^ 3T. (36) 

As Kn gets large, Tan can get small even in the weak coupling limit. Thus, 
large gradients, which are the hallmark of the Neck zone, increase the impor- 
tance of anomalous relaxation over color diffusion and collisional relaxation. 
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Short relaxation times arise not only in the strong coupling but also in the 
weak coupling but classical field limits. 

The full source (t^{X) was computed analytically by Neufeld in Ref. ¡60] 
and it reads 

a'^(X) = ((7°,v(tO -a^). (37) 
The components of the source term are reproduced below for convenience 



(T°(í,x) = (¿(í,x)7i;^ 



Xi 



{xf + xj_) 



Xl + 



7 vx. 



x\ + xy^ 



(38) 



cj4í,x) = (x-ví) ^ 



87r {xl + xlY 



4 4 
V X , 



+ {xy^ + xl) {2x1 + {v'^ + 2) xl/-f^) 2vxi 



{xy^ + xif ^ Uy^ + 



(39) 



where 7 = 1/Vl — «s = ff^/ (47r) = l/vr, = and C = 4/3 for 
a quark. These equations were obtained in the hmit where the dielectric 
functions that describe the medium's response to the color fields created by 
the heavy quark were set to unity. The eíFects from médium screening on 

were studied in detall in Ref. [6OJ . In our numerical calculations we used 
x±max = ^/ttT'D as an infrared cutoff while the mínimum lattice spacing 
naturally provided an ultraviolet cutoff. The background temperature was 
set to To = 0.2 GeV. 

This result holds for the weak coupling, long relaxation time 1/r* = 
e ^ O"*" limit. In this limit the color conductivity is large and Joule heat- 
ing efficiently converts field energy and momentum into plasma heating and 
collective flow. The most interesting aspect of this for the present applica- 
tions is the effect of highly inhomogeneous anisotropic Lorentz contracted 
Coulomb field of a uniformly moving heavy quark and color Casimir 4/3 on 
the collective fiow pattern imprinted on the plasma. 



4-2. Numerical Results for the Freeze-out 

We now turn to the observable consequences of the pQCD chromo-fiuid 
flow. The plot in Fig. [12] was obtained using the source of section 14.11 
integrated into the hydrodynamical code used in section [2l The associated 
heavy quark jet is created in the beginning of the hydro evolution í = O 
at Xl = —4.5 fm and the freeze-out is done when the it r caches the origin 
of the coordinates, independently of the heavy quark's velocity (note that 
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x = z-vt (^/KTJ 

Fig. 12. (Color online) The relative local temperature perturbation AT/Tq = 
T(xi, x±)/Tq — 1 and flow velocity profile due to a heavy supersonic quark jet 
moving with v — 0.9 (speed of sound Cs = 1/VS)- The results were obtained us- 
ing perfect fluid (3+l)D hydrodynamics in the presence of the pQCD source term 
computed by Neufeld in [ÜD]. The panel shows the Mach wake (see purple dashed 
line) and trailing shear column in the far zone as well as the Neck región (red) near 
the jet. The heavy quark is at the origin of the coordinates. The arrows show the 
direction and magnitude of the flow. The numbers in the plot label the contours of 
constant AT/Tq. Note that non-Mach flow induced by Joule heating is generated 
near the jet. 



we do not include the trigger jet in our analysis). This provides a very 
rough description of the case in which a very energetic heavy quark punches 
through the médium. 

Comparing Fig. [T2]with Fig.[Ü]it is apparent that the qualitative features 
displayed by both of these systems are very similar. In particular, both 
Systems exhibit the Mach-like behavior when the distance from the source 
~ 5/T. The transverse ñow in the Neck región, however, is considerably 
smaller in pQCD in comparison to the AdS result found in section 13.31 We 
then expect, based on the study in section 13.31 that the pQCD produced 
Mach cones are not observable. 

We present normalized CF azimuthal distribution in the form of Eq. 
[5] in Fig. [13l We computed /((/>) for v = 0.58,0.75,0.90 in a static uni- 
form background. The results for the azimuthal angular correlations for 
Pt = SvrTo = 5 GeV, and y = O are shown in Fig. [13l The pQCD angu- 
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Fig. 13. (Color online) Normalized and background subtracted azimuthal away 
side jet associatcd correlation after Cooper-Frye freeze-out (Eq. [S]) for pQCD. 
Here /(0) = cfiN/dydp^díp is evaluated at p± — 12.5roj and y = 0. The black line 
is for V = 0.58, the blue line for v — 0.75, and for the red line v = 0.9. Compare to 
the AdS/CFT rcsult in Fig.[llJ 



lar distribution shows only a sharp peak at vr for all velocities while the 
AdS/CFT distribution displays the double-peak structure for all velocities 
shown in Fig. [TTl Note, once again, that the peaks in the AdS/CFT cor- 
relation functions do not obey Mach's law since they do not come from the 
far-away linearized región. The smaller transverse flow of the Neck produced 
by the pQCD source term ensures that the Mach contribution is washed out 
by the Cooper-Frye thermal smearing. 



5. Conclusions 

In conclusión, we have shown that Mach cones in the fluid phase of the 
system are a generic prediction of a médium characterized by low viscosity 
and a high opacity to hard probes. However, the signal seen in the ex- 
periment is likely to be considerably different from the naive expectation 
both because of the medium-induced thermal smearing at freeze-out and 
the presence of additional structures such as the diffusion wake. 
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We have illustrated these issues using examples taken from AdS/CFT 
and pQCD and shown that the critical requirement for a "Mach-cone like 
signal" is in fact the presence of a strong transverse ñow in the non-linear 
unthermalized Neck región cióse to the source since this is most likely to 
survive freeze-out and yield a cone-hke signal. However, the formation of 
strong transverse ñow in the Neck región is far from assured. For instance, 
it is prominent in the AdS/CFT solution and missing in pQCD. 

Our analysis is far from complete. The most obvious physical aspects it 
misses are the underlying transverse ñow (probably significant at timescales 
relevant to jet absorption), coalescence dynamics (certainly signiñcant at 
momenta comparable to the away-side correlation used by experiment), and 
the dynamics of the freeze-out hypersurface. Efforts in this direction are on- 
going [63] and more studies are needed before a conclusive link between the 
experimental result and the theory can be made. Whille the Mach cone 
interpretations [26] of conical correlations are intriguing and have focused 
attention on modérate pT correlation observables, the observed conical cor- 
relations will only fulfill its promise as a unique probé of the sQGP dynamics 
and its coupling to hard probes when the many remaining open problems 
and caveats discussed above are resolved. This effort is certainly warranted 
given the potential power of such correlation observables to resolve dynam- 
ical features inaccessible with well studied and measured inclusive and bulk 
ñow observables. 
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